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Abstract 

Some special features of magnetically deformed multi-electron atoms in the 
crustal region of strongly magnetized neutron stars (magnetars) are studied 
with the relativistic version of semi-classical Thomas- Fermi model. An analyt- 
ical solution of Poisson's equation with cylindrical symmetry is obtained. The 
variation of atomic dimension with magnetic field strength is also studied. An 
expression for atomic quadrupole moment, appears because of magnetic de- 
formation is derived. The effect of atomic quadrupole-quadrupole interaction 
in the equation of state of crustal matter of magnetar is briefly discussed. 
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From the observational evidence of some compact stellar objects with extremely high 
surface magnetic field 0, the study of the effect of ultra strong magnetic field on the geo- 
metrical structure (electronic orbits) of atoms in the crustal region have got a new dimension. 
The possibility of such ultra-high magnetic field in these stellar objects has given re-birth 
to this "the old subject" of quantum mechanics. These strongly magnetized objects are 
known as magnetars fl2|. The measured value of magnetic field strength at the magnetar 
surface can be as high as > 10 15 G. These objects are also supposed to be the sources of soft 
gamma repeaters and anomalous X-ray pulsars. Soft gamma repeaters are a class of com- 
pact stellar objects which periodically emit bursts of gamma rays, in addition to persistent 
X-rays. These two classes of objects appear to be very similar (aside from the presence of 
bursts in SGR's). The measured period lie in the range 5 — 12 sees. The dipole fields at 
the surface inferred from the spin down rate are 10 14_15 G, which are much larger than the 
surface field strength of known strongly magnetized conventional neutron stars (~ 10 12 G). 
It has long been suggested that strong magnetic energy may be the ultimate source of both 
the bursts and the persistent radiation [Q. From the observational data it is believed that 
such objects (magnetars) are strongly magnetized young neutron stars. If the strength of 
surface magnetic field is > 10 15 G, then it would be interesting to investigate the effect of 
such strong magnetic field on the properties of crustal matter. Some of the properties of 
such low density matter in presence of strong magnetic field have already been studied ||. 
The other interesting effect of ultra-strong magnetic field on the crustal matter, which has 
not been explored in greater detail, which we are going to discuss in this article is the possi- 
bility of magnetic deformation of atoms at the lattice points in the crustal region. Although 
the physical problem related to the atoms placed in an external magnetic field is one of 
the oldest subject of quantum mechanics, this particular exotic phenomena - the magnetic 
deformation of atoms, heavier than hydrogen, have not been studied. The effect of strong 
magnetic field on the stability of hydrogen atoms and on the multi-proton bound states have 
been investigated in detail 0-0. In a recent work, the magnetic deformation of hydrogen 
atom and the appearance of its quadrupole moment have been studied [f|. In the present 
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report, we shall consider the effect of strong magnetic field on the geometrical structure of 
multi-electron atoms (much heavier than hydrogen atom). We do believe that such physical 
situation can not be ruled out at the crustal region of magnetars. We shall also investigate 
the effect of strong magnetic field on the quadrupole moment of deformed atoms at the 
crustal lattice points. It is believed that the shape of an atom is almost spherical at B < B 
and a cylindrical at B > Bq, where Bq = m^e 3 = 2.35 x 10 9 G, m e and e are respectively the 
mass and electric charge of the electron (in our calculation we have assumed h = c = 1). 
The critical parameter 7 = B/Bq has significance in the non relativistic theory. In our case 
B > 4.4 x 10 13 G which is the relativistic region and we consider the parameter ( = B/B c , 
where B c = 4.4 x 10 13 G, the typical field strength at which the Landau levels of electrons are 
populated. Since fully relativistic calculation of a many body quantum mechanical system of 
finite size in presence of ultra-strong magnetic field is a very difficult task, even numerically, 
we are considering for the multi-electron atoms in presence of strong magnetic field B > B c , 
a simple semi-classical approach- the relativistic version of Thomas- Fermi model in presence 
of strong magnetic fields. 

We shall now consider the basic formalism of relativistic Thomas- Fermi model in presence 
of strong magnetic field. The electrostatic potential (j) satisfies the equation Q], 



Where n e is the electron density, Z is the atomic number and r n is nuclear radius. Since 
r n rsj fm <C atomic radius ~ A , the nuclear contribution on the right hand side can therefore 
be neglected for r > r n . Now assuming that only the zeroth Landau level for the electrons 
is populated, the chemical potential of electron is given by 



V 2 = Airn e e — AnZe5(r — r n ) 



(1) 




(2) 



This gives, 




(3) 
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the electron Fermi momentum. Now in presence of a strong quantizing magnetic field of 
strength B, the electron number density is given by, 



7T Z 7T 



In the high density region, we can neglect m e , then 



and eqn.(l) becomes 



n e ~^(/i + e0) (5) 



V 2 <j ) =—( f i + e<f ) ) (6) 
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Substituting jj, + ecp — ip, we have in the cylindrical coordinate (r, 9, 0) 

iprr + - A + + ipzz = ail> (7) 

Where, ip x = dtp/dx and 

a = (8) 

7T 

Assuming cylindrical symmetry, becomes independent of 9. Then with the separable form 
of solution 

^(r-z) = R(r)Z(z) (9) 

we have 

rP 7 

" - s 2 Z = (10) 



and, 



dr 2 r dr 

Where, p 2 = s 2 = a/2. The solution of the first equation is 



dz 2 



d 2 R 1 dR 9 , . 

+ -3 p 2 R = (11) 



= constant x exp(— sz) (12) 
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Replacing r by r = pr, we have from the second equation 

d 2 R 1 dR 



dr 2 ~^ r dr ^ ^ 

The solution of this equation is given by, 

R(r) = constant x K (f) (14) 

where K (f) is the second kind modified Bessel function of order zero. Therefore, 

ip(r,z) = C x exp(-sz)K (r) = C x exp(-z)K Q (r) (15) 

To obtain the normalization constant C, we use the relation, 

lim ip = lim e<pl = Ze 2 (16) 

Where / = (r 2 + z 2 ) 1 ^ 2 , some arbitrary distance from the centre of the nucleus to the surface 
of the cylinder. Since r n r atom , we can approximate the normalization constant by 

Ze 2 

exp(f B ) (17) 

r n K (r n ) 

which is ~ 10 11 for iron atoms and is insensitive with the variation of magnetic field strength. 
Since r = f/p and z — z/s, we can conclude that both the longitudinal and lateral dimensions 
of the deformed atom decrease with the increase in magnetic field strength. To obtain the 
surface of a deformed atom we use the relation, 

which means the electric field vanishes at the surface. Hence we get, 

exp(-z)[fK (f) + zK^f)} = (19) 

Solving this equation we obtain r as a function of z. We have noticed that the atoms become 
ellipsoidal with cylindrical symmetry in presence of strong magnetic field. The variation of 
r with z is given by the equation 



z 2 r 2 



a 



+ ^ = 1 (2°) 



The length of both the major and minor axes decrease with the increase of field strength. 
The variations are given by the following power laws. 



a = 0.53C" - 5 (21) 



b = 0.436C U - 5U1 (22) 

From these power law relations it is very easy to see that the decrease of b with B is slower 
than that of a. This is also obvious from the asymptotic nature of Kq(v) which is given by 



W ~ « (23) 

for large f. Which means the deformed atom maintains its ellipsoidal shape even for very 
high value of magnetic field strength. The quadrupole moment for the deformed atom with 
ellipsoidal shape of the type given by eqn.(20) may be written as [p|, |lO| , [TT 



Qi j = ^p{2x i x j -l 2 S ij ) (24) 

where the sum goes over all the charges. Now for the charge distribution with cylindrical 
symmetry 

Qzr = Qrz = (25) 

and 



Q zz = -Q„ = 7T / p(z, r){z 2 - r 2 )r 2 dz (26) 



where, 



eB 

P(r,z) ~ —^(r,z) (27) 



7T 2 
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Expressing r in terms of z using eqn.(20), we have obtained Q zz as a function of (. The 
variation of Q zz with ( is given by the power law, 

Q zz = Cx 0.029C 1 ' 5 (28) 

for 7 >> 1, C is the same normalization constant as mentioned before. However, for a very 
broad range of (, the dependence is completely different, particularly for low (, In fig.(l) 
we have shown the change of Q zz as a function of ( for a very wide range. The variation 
is almost identical with the recent result on deformed hydrogen atom reported by Potekhin 
§• 

We shall now consider a deformed atom in presence of external electric field produced 
by the quadrupole moment of other deformed atoms. Let us assume that there are two such 
systems, each having total charge zero and vanishing dipole moment. Then we have the 
total potential energy of the system 

C/ = £><Mr!) (29) 

i 

where fl is the radius vector of the charge e« with its origin anywhere within the system. 
Since both the total charge and dipole moment of the system are zero, we have the total 
potential energy due to quadrupole field 

where (p^ is quadrupole field and Q a p is the quadrupole moment. The quadrupole field is 
given by, 

4* - (3D 

which is the quadrupole field produced by the quadrupole moment Q a s of a deformed atom 
at the point with radius vector R Q . Combining these two equations we have after some 
straight forward algebra, 



This is the potential energy of the deformed atom in presence of the quadrupole potential (f)^ 
produced by another deformed atom, i.e., it is the extra two body potential arises because 
of quadrupole-quadrupole interaction. 
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FIGURES 
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FIG. 1. Variation of quadrupole moment (in A 2 ) with the parameter ( = B/B c , the scaling 
parameter is the normalization constant C ~ 10 11 . 
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